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Abstract 

A Woods-Saxon equivalent to a double folding potential in the surface 
region is obtained for the heavy-ion scattering potential. The Woods- 
Saxon potential has fixed geometry and was used as a bare potential 
in the analysis of elastic scattering angular distributions of several sta¬ 
ble systems. A new analytical formula for the position and height of the 
Coulomb barrier is presented, which reproduces the results obtained using 
double folding potentials. This simple formula has been applied to esti¬ 
mate the fusion cross section above the Coulomb barrier. A comparison 
with experimental data is presented. 


1 Introduction 

The strong nuclear force is responsible for keeping the nucleons together inside 
the nucleus and is still not fully understood. The attractive force between the 
nucleons is the residuum of the interaction between quarks and gluons confined 
inside the nucleons and the connection between the fundamental interaction and 
the nucleon-nucleon force is still an open problem. Elastic scattering between 
nuclei provides information of the nuclear interaction however, the cross sections 
are affected by the couplings between the elastic scattering and all other possible 
reaction channels. The potential obtained from the analysis of elastic scattering 
angular distributions is the sum of a bare potential, which is real in principle, and 
a polarization term, which is complex and contains the effects of all couplings. 
To obtain information of the bare potential one should find a situation where 
the elastic is the only open channel, however, it is very difficult to find such 
experimental situation. Even at energies around the Coulomb barrier, where the 
reactions channels are closing, there is still the contribution of the fusion process, 
which makes the interacting potential complex. If we go down to even lower 
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energies, the effect of the short range nuclear potential becomes smaller and 
smaller as the long range Coulomb potential dominates, making the scattering 
pure Rutherford. 

The bare potential can be in principle be defined as the result of the double 
folding of the nucleon-nucleon interactions and the projectile and target nuclear 
densities [TJ |21 131 SI- The nucleon-nucleon interactions can be obtained from 
more fundamental theories. Double folding potentials have been used as the 
bare potential to analyse experimental data [S] and the imaginary part of the 
interaction is normally parameterized and freely searched to best fit the angular 
distributions. One of the most widely used parameterizations for the nuclear 
potential is the well known Woods-Saxon (WS) shape [5], with three parameters 
that are adjusted to reproduce the data. 

More recently, the Sao Paulo optical potential (PSP) has been developed, 
where the real part is taken as a double folding potential and the imaginary 
part has the same geometry of the real part, with an additional fixed normal¬ 
ization factor. An energy dependence term is included to account for non-local 
corrections due to the Pauli principle [Z1[B]. The Sao Paulo potential has no free 
parameters and has been succesfully applied to a large number of experimental 
angular distributions from low to intermediate energies. 

Despite the success of the Sao Paulo potential in analysing elastic scattering 
data, it would be interesting to investigate the relation between the Woods- 
Saxon shape and the double folding potential. Most of the optical model and 
reaction programs use the WS parameterization whereas double folding poten¬ 
tials have to be entered externally from numerical files. The equivalence between 
WS and double folding potentials is not straightforward and there is no WS that 
could reproduce the double folding shape in the whole radial range. However, 
heavy ion scattering at low energies is frequently sensitive only to the tail of 
the nuclear potential and not very much to the potential in the interior region, 
where strong absorption usually takes place. 

In the present work we determine the parameters of a WS potential that 
reproduce the real part of the double folding Sao Paulo potential in the surface 
region. We determine the geometry of this potential and apply it to experimental 
elastic scattering angular distributions of several systems. 

Based on the real part of this potential we found a simple analytical formula 
to obtain the position and the height of the Coulomb barrier which reproduces 
quite precisely the Coulomb barriers obtained from the Sao Paulo potential. 
This formulation is applied to estimate the fusion cross section in the region 
above the Coulomb barrier. 


2 The Woods-Saxon potential 

The Woods-Saxon shape is given by: 


f{r) 


1 

l + exp(^) 


( 1 ) 
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Figure 1: Radial dependence of the real part of the nuclear potential for several 
systems. 

where R = and a is the difuseness. The optical potential is 

written as: 


Vnuci{r) = -Vofi{r) - iWof2ir) (2) 

where Vq and Wq are the real and imaginary strengths of the optical potential 
respectively and fi{r) are Woods-Saxon form factors that may have different 
values of radius and diffuseness parameters. 

We developed a simple computer program to perform an automatic search 
on the 3 parameters Vq, ro, a oi Eq. [T]to reproduce the tail of the real part of 
the Sao Paulo potential for several systems. 

The results are shown in Fig. [T] where the dashed lines represent the real 
part of Sao Paulo potential (PSP) calculated for several systems. The solid lines 
are the resulting Woods-Saxon potentials that best fit the PSP in the region 
r > R. We found that, for all systems analysed, the Woods-Saxon that best 
reproduces the tail of the PSP potential has a diffuseness a ss 0.65 fm, and 
ro ~ 1.3 fm for strengths ranging between Vq = 10 — 20 MeV. In general, double 
folding potentials are strongly attractive with strengths of hundreds of MeV’s. 

2.1 Ambiguities 

Ambiguities in the optical potential have been the subject of many studies 
since the beginnings of nuclear physics [5]. The fact that the elastic scattering 
angular distributions in the strong absorption regime are sensitive to a small 
region in the surface of the nuclear potential, has long been recognized. If the 
scattering at low energies is sensitive only to the tail of the nuclear potential, 
for a Woods-Saxon shape one immediately gets that, for (r — R) » a, we 
obtain Vnuciir) = —Vq exp(i?/a) exp(—r/a). Thus, for a given difuseness a, any 
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combination of Vq and R that leaves Vbexp(i?/a) unchanged will provide the 
same potential in the surface region. In this sense the parameters Vq and tq 
proposed in the present paper are only one family of possible potentials. 


3 The Coulomb barrier 


The height and position of the Coulomb barrier are very important parameters 
in the collision of two heavy ions. They basically determine the total reaction 
cross section at energies above the Coulomb barrier and can be obtained if the 
real part of the nuclear potential in the surface region is known. 

The condition: 

-r[^nucl{r) + Vcoul{r)]r=RB = 0 (3) 

dr 

determines the position and the height of the Coulomb barrier is given by 
Vb = VnuciiRs) + Vcoui{Rb)- As usually Vnud is unknown, we may assume an 
approximate radius for the Coulomb barrier radius as: 


RB=R = ro{A\/^ (4) 

and use a simplified formula for the height of the Coulomb barrier: 


Vb = 


Z\Z2^ 

Rb 


(5) 


In general Eqs. U and 0 do not yield very good results in comparison to 
Coulomb barriers obtained from realistic DF potentials. Eq. |3] provides a posi¬ 
tion for the Coulomb barrier which is, in most cases, larger than the geometrical 
radius from Eql4] provides similar values of i?_B. Indeed, we see that, if we take a 
Woods-Saxon form for the nuclear potential Vnuci{r) = —Vo/[l+exp[(r —i?)/a]] 
and Vcoui{t) = Z\Z 2 ^jr in Eq. [3] one shows that, in the approximation 
exp((i?B — R)la) » I, the Coulomb barrier radius can be written as: 


Rb = i? -b a In 


[A K)] 

— X 

. ® Vb. 


( 6 ) 


This equation may not be exact but it displays the main physics of the 
relation between R and Rb- For a square potential (a = 0), we get Rb = R 
and there is no correction to i?. As the difuseness of the potential increases, the 
correction term in right hand side of Eq. [5] increases. Also, for larger Coulomb 
barriers, the correction decreases. 

Taking R — -b fm, a = 0.65 fm, Vq = 15 MeV and Vb given 

by Eq. [S]one obtains: 


Rb = R + 0.65 ln[a;] 


where x = 27.1 x 




1 / 3-,2 
2 ) 


Zi Z2 


is a positive dimensionless parameter: 


(7) 
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Then we get for Vb' 


(8) 


Rb a: + 1 

Equations [7] and [5] depend only on the masses and charges of the nuclei and 
provide the Coulomb barrier position and height in very good agreement with 
those obtained from numerical calculations using the Sao Paulo Potential. This 
is shown in Table [T] The discrepancies are in the order of a few percent or 
smaller than that for heavier systems. 


3.1 The curvature of the Coulomb barrier 


One could go a step further and obtain the curvature of the Coulomb barrier 
based on the above potential. The region around the top of the Coulomb barrier 
can be approximated by an inverted harmonic oscillator potential of height Vb 
and frequency w. The frequency is related to Vb by: 


hw = h^J {\(FV (r) / dr'^\)r=RB / ^^ 
where V{r) = Vnuciir) + Vcoui{r), is the reduced mass and: 

Vo X{X - 1) 2^1^262 


{d^V{r)/dr^)r=Rs = - 


{x + 1 )^ 


R3 


(9) 


( 10 ) 


Substituting Vq = 15 MeV and a = 0.65 fm in the above formula, one can 
estimate the curvature of the Coulomb barrier. 


Table 1: Comparison between formulas [7] and H] (WS) and the results for 
Coulomb barrier radius and height from numerical calculations using the double 
folding Sao Paulo potential (PSP). 


(fm) VB(MeV) 


System 

X 

WS 

PSP 

WS 

PSP 

“He+^Be 

51.45 

7.62 

8.00 

1.22 

1.32 

12C+160 

13.06 

7.92 

8.15 

7.65 

7.78 

160+24Mg 

8.24 

8.39 

8.55 

14.84 

14.85 

16o+58Ni 

4.94 

9.34 

9.40 

31.99 

31.68 

9Be+208pb 

5.29 

11.48 

11.50 

38.72 

38.55 

16o + 208pb 

2.94 

11.68 

11.65 

77.07 

75.90 


4 Analysis of experimental data 

In the next sections we use this potential to analyse a number of experimental 
angular distributions. We fix the geometry of the real and imaginary parts to 
ro = 1.3 fm and a = 0.65 fm and allow Vb and Wo vary to best fit the angular 
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distributions. The idea here is not to obtain excelent fits with such a simple 
potential but to show that it is possible to reproduce the general features of the 
angular distributions using this fixed geometry potential. 

All the Optical Model calculations have been performed using the program 
SFRESCO, the automatic search version of the FRESCO program [TU]. The 
results are presented in the next subsections. 

4.1 ®Be+^^Al system 

Nine elastic scattering ®Be+^A.l angular distributions measured by P. R. S. 
Gomes et al. m have been analysed in the range from 12 MeV to 40 MeV 
in the laboratory system. Equation [5] gives Vb = 10.65 MeV for the Coulomb 
barrier energy in the laboratory system. The geometry of the real and imaginary 
potentials was fixed at Tq = 1.3 fm and a = 0.65 fm and the depths Vq and Wq 
were varied to best ht the data. The Coulomb radius parameter was fixed at 
roc = 1.3 fm. The results are shown in Figure [2] and the fitted parameters are 
presented in Table [2] together with the errors and the best reduced chi-square 
values. The errors have been estimated by the program SFRESCO using the 
gradient method used in the search procedure. 


Table 2: Parameters of WS potential for the system ®Be-|-^Al. 


Eiab (MeV) 

Vo (MeV) 

Wo (MeV) 


12.0 

18.15± 2.10 

15.30± 7.68 

0.36 

14.0 

15.22± 1.46 

10.08± 2.90 

0.47 

18.0 

17.87± 0.21 

3.37± 0.17 

4.35 

22.0 

16.45± 2.72 

16.50± 4.05 

0.69 

25.0 

12.22± 0.30 

13.64± 0.01 

1.83 

28.0 

17.39± 0.75 

12.17± 0.52 

3.92 

32.0 

15.90± 0.62 

11.96± 0.52 

6.24 

33.0 

17.43± 0.48 

12.12± 0.35 

13.36 

40.0 

8.29± 1.33 

10.62± 1.85 

2.70 


4.2 system 

Eleven angular distributions from 35 MeV to 48 MeV have been analysed [T^ . 
Vj^b _ 4 Q gQ MeV for this system. The results are presented in Figure |3] and 
in Table [3] The small values and large errors of Vq shown in Table |3] for Eiab = 
35 — 37 MeV show that the potential is not determined at these energies below 
the Coulomb barrier. 
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Figure 2: Cross section for the ®Be+^^Al system. The energies values are shown 
in tabled] The experimental data is from Ref. |13) . 
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Figure 3: Angular distributions for the ^®0+®®Ni system. 
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Table 3: Parameters of WS potential for the system ^®0+®®Ni. 


Slab (MeV) 

Vo (MeV) 

Wo (MeV) 


35.0 

O.lOi 4.35 

3.84± 0.41 

1.27 

35.5 

O.lOi 1.51 

3.35± 0.27 

1.51 

36.0 

O.lOi 3.67 

3.54± 0.35 

0.91 

36.5 

2.56± 4.43 

2.63± 0.80 

2.01 

37.0 

2.00± 0.06 

3.48± 0.60 

1.98 

38.0 

15.95± 0.79 

l.OOi 0.19 

0.66 

40.0 

17.59± 0.002 

1.18± 0.09 

1.06 

44.0 

11.79± 0.0006 

7.61± 0.004 

9.48 

46.0 

11.05± 0.007 

9.68± 0.009 

13.17 

48.0 

11.42± 0.0002 

ll.lli 0.0005 

9.03 


4.3 ®Be+®'^Zn system 

Six angular distributions from 17 MeV to 28 MeV in the laboratory system 
have been analysed [14]. The Coulomb barrier is at 19.4 MeV in the laboratory 
system (Eq. [8|) and the results are presented in Figure |4| and the resulting 
parameters are in table S) 


Table 4: Parameters of WS potential for the system ®Be+®^Zn. 


Eiab (MeV) 

Vo (MeV) 

Wo (MeV) 

17.0 

0.000 ± 

4.90 

36.21± 

3.17 

19.0 

12.09± 

1.43 

19.17± 

1.75 

21.0 

12.06± 

0.29 

18.88± 

0.61 

23.0 

10.61± 

0.16 

14.51± 

0.40 

26.0 

10.12± 

0.18 

15.98± 

0.79 

28.0 

13.01± 

0.25 

12.47± 

0.49 


4.4 ®Be+^®Y system 

Seven angular distributions have been analysed for the ®Be+®®Y system. The 
laboratory energies range from 19 MeV to 33 MeV and = 23.25 MeV. 
The results are presented in Table [5| The values of the imaginary part of the 
potential drop down to energies lower than the Coulomb barrier, except for the 
lowest energy. 
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Figure 4: Cross section for the ®Be+®^Zn system. The energy values are shown 
in table m The experimental data is from Ref. |14) . 
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Figure 5: Cross section for the ®Be+®®Y system. The energy values are shown 
in table [5] The experimental data are from Ref. [T5] . 
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Table 5: Parameters of WS potential for the system ®Be+®®Y. 


Eiab (MeV) 

Vo (MeV) 

Wo (MeV) 


19.0 

13.26± 0.40 

17.37± 0.40 

2.70 

21.0 

20.51± 0.92 

1.36± 0.87 

1.87 

23.0 

10.07± 0.68 

7.48± 0.71 

1.51 

25.0 

8.50± 0.01 

7.93± 0.01 

2.24 

27.0 

2.52± 0.34 

29.12± 0.72 

17.46 

29.0 

8.63± 0.16 

11.59± 0.33 

11.81 

33.0 

14.42± 0.10 

15.11± 0.13 

11.74 


4.5 ®Be+^^^Sm system 

Ten angular distributions for the ®Be+^"'^'^Sm system have been analysed [ISIIS]. 
The energies range from below to above the Coulomb barrier at = 33.24 
MeV. See Figure [Hand Tabled 


Table 6: Parameters of WS potential for the system ®Be+^^^Sm. 


Eiab (MeV) 

Vo (MeV) 

Wo (MeV) 


30.0 

11.99± 2.55 

15.28± 1.09 

1.16 

31.5 

13.94± 0.81 

12.76± 0.57 

1.32 

33.0 

0.81± 4.24 

16.76± 2.47 

0.98 

34.0 

7.44± 1.59 

13.03± 1.62 

0.73 

35.0 

9.97± 0.83 

13.41± 1.31 

0.92 

37.0 

7.42± 0.58 

17.00± 1.11 

0.96 

39.0 

11.86± 0.42 

13.41± 0.78 

2.58 

41.0 

11.63± 0.45 

15.81± 1.04 

1.87 

44.0 

13.66± 0.18 

15.67± 0.41 

1.26 

48.0 

13.83± 0.16 

16.86± 0.35 

14.34 


4.6 ®Be+^°®Pb system 

Fourteen angular distributions have been analysed in an energy around and 
below the Coulomb barrier [T7]. = 40.40 MeV for this system. 


5 Fusion and total reaction cross section 

The total reaction cross section is an important information that can be obtained 
from the elastic scattering. In addition, at energies around the Coulomb barrier, 
in many cases, fusion exhausts most of the total reaction cross section and can 


12 









Figure 6: Cross section for the ®Be+^^'^Sm system. The energies values are 
shown in table [HI The experimental data are from Ref. [161 [T5] . 
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Figure 7: Cross section for the ®Be+^°®Pb system. The parameters are shown 
in tableland the experimental data are from Ref. |17) . 


be estimated by barrier penetration calculations. The well known Wong formula 
[18] for fusion has been applied with success to provide estimations of the fusion 
cross section. 


^wong — 


R% hw 
TR 


ln{l + exp[ 27 r(i? — VB)/hw\} 


( 11 ) 


This formula depends on 3 parameters, the Coulomb barrier position, height 
and its curvature {hw), the same parameters that have been determined on 
Sec.III. For energies above the Coulomb barrier the Wong formula reduces to a 
simpler one which depends only on two parameters, the position and height of 
the Coulomb barrier. 

a = ^Rl{l-^) ( 12 ) 
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Figure 8: Comparison between the total reaction, fusion cross section with 
predictions of formula [TT] and [12] for the ®Be+^^Al system with Rb = 8.29 
fm, = 7.99 MeV and hw = 3.25 MeV. The fusion cross sections were 

taken from ref. m and the reaction cross section are from the optical model 
calculations from section KT\ 



Figure 9: Comparison between the total fusion cross section with predictions of 
formula [TT] and IT^ for the ®Be+®^Zn system with Rb = 9.28 fm, = 17.0 

MeV and Hw = 3.44 MeV. The fusion cross sections are taken from ref. El¬ 


is 









Figure 10: Comparison between the total fusion cross section with predictions 
of formula ITT] and IT^ for the ^"^N+^^Co and ^®0+®®Co systems. The fusion cross 
sections are taken from ref. |16) . 



Figure 11: Comparison between the total fusion cross section with predictions 
of formula m and HD for the ®Be+^'^^Sm system | 


Rb = 10.67 fm, vq 


yicm) ^ 


31.29 MeV and hw = 3.54 MeV. The fusion cross sections are taken from ref. 
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Table 7: Parameters of WS potential for the system ®Be+^°®Pb. 


Eiab (MeV) 

Vg (MeV) 

Wo (MeV) 


37.0 

0.10 ± 2.77 

29.38± 0.56 

6.34 

37.8 

13.95± 0.84 

21.58± 0.48 

2.12 

38.0 

15.98± 0.64 

19.98± 0.39 

3.37 

38.2 

17.52± 0.54 

18.69± 0.36 

3.91 

38.5 

18.39± 0.37 

18.53± 0.28 

3.59 

38.7 

18.01± 0.53 

17.27± 0.42 

3.34 

39.0 

19.08± 0.25 

17.18± 0.22 

7.60 

39.5 

18.20± 0.20 

17.79± 0.20 

5.02 

40.0 

16.68± 0.18 

19.06± 0.20 

7.07 

41.0 

14.84± 0.07 

19.76± 0.01 

9.17 

42.0 

13.01± 0.10 

20.69± 0.14 

18.00 

44.0 

11.94± 0.19 

20.89± 0.28 

8.01 

47.2 

14.50± 0.07 

19.84± 0.18 

13.60 

50.0 

15.15± 0.11 

22.81± 0.24 

28.47 


We applied formulas llll and ll2l using the parameter calculated from formulas 
[7l [3 [9] and uni for the ®Be+^^Al, ®Be+®^Zn and ^®0+^°®Pb systems. The results 
are shown in figures El 13 [13 El and [13 

We see that formulas [TT[ (solid) and [T3 (dashed) give the same result as 
the energy overcomes the Coulomb barrier. The agreement between the exper¬ 
imental fusion cross section and calculation is reasonable for energies above the 
Coulomb barrier. For energies below the barrier the calculations predict a cross 
section much smaller than the experimental one as can be seen in Figure 1101 
This is expected since it is well known that for energies below the Coulomb 
barrier the fusion is strongly affected by coupled channels effects, which are 
obviously not taken into account by the simple formulation presented here. 


6 Conclusions 

The real part of a Woods-Saxon potential that fits a double folding potential in 
the surface region (r > R) has been obtained. It was found that the tail of the 
double folding potential can be very well reproduced in all cases analysed here 
using a Woods-Saxon potential with a fixed geometry rg = 1.3 fm and a = 0.65 
fm and depths varying between 10 — 20 MeV. There is a continuum ambiguity 
between rg and Vg. 

A simple analytical formula has been derived using the real potential with 
depth Vg = 15 MeV, which provides the position and height of the Coulomb 
barrier in very good agreement with the double folding potential predictions. 
It is shown that the Coulomb barrier position and height depend on a single 
dimensionless parameter x, which can be easily calculated as a function of the 
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Figure 12: Comparison between ^®0+^°®Pb fusion cross sections m and the 
calculations using the Wong formula and Rb = 11.68 fm, = 77.07 MeV 

and hw = 2.45 MeV. 


masses and charges of the colliding nuclei. 

An optical model analysis has been performed for several systems using such 
potential, fixing the geometry for the real and imaginary parts, and adjusting the 
depths V and W to fit the angular distributions. It is shown that the potential 
proposed here provides reasonable fits of the scattering angular distributions for 
several stable systems at several energies above and below the Coulomb barrier. 
A strong variation of Vb and W is observed at energies around the Coulomb 
barrier, as expected, due to the closing of the reaction channels as the energy 
goes down below the Coulomb barrier. 

A criticism could be done to the optical model analysis presented here since, 
the obtained potentials, are not anymore strictly equivalent to the double fold¬ 
ing, not even in the suface region, because of the free variation of the real and 
imaginary depths. This is true and, it must be indeed just like that since the 
optical potentials that reproduce the data are not anymore the bare potential 
but the total optical potential, which includes all the polarization effects from 
the couplings with other reaction channels. 

Total reaction and fusion cross sections have been calculated using the ana¬ 
lytical formula derived here and the result is compared with experimental fusion 
cross sections. It is shown that, above the Coulomb barrier, the analytical for¬ 
mula provide a good approximation for the fusion cross sections and the calcula¬ 
tions can be done with a simple pocket calculator without the need of numerical 
computations. 
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